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Abstract 

In this paper, we present standard Johnson UBV photometry of the echpsing 
Q-[ binary BD +36 3317 which is known as a member of Delta Lyrae (Stephen- 

son 1) cluster. We determined colors and brightness of the system, calculated 
E{B — V) color excess. We discovered that the system shows total eclipse 
in secondary minimum. Using this advantage, we found that the primary 
component of the system has 59 — AO spectral type. Although there is no 
pubhshed orbital solution, we tried to estimate the physical properties of the 
system from simultaneous analysis of UBV light curves with 2003 version of 
Wilson - Devinney code. Then we considered photometric solution results to- 
gether with evolutionary models and estimated the masses of the components 
as Ml = 2.5 Mq and M2 = 1.6 Mq. Those estimations gave the distance of 
the system as 353 pc. Considering the uncertainties in distance estimation, 
resulting distance is in agreement with the distance of Delta Lyrae cluster. 
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(BD+36 3317) 



1. Introduction 



BD+36 3317 (SAO 67556, asooo = IS'' 54"* 22^ (Jaooo = 36° 51' 07", V = 
8"*. 77) is an Algol type eclipsing binary which i s in the same area with D elta 



Lyrae cluster (02000 = 1 



2000 



36° 49', 



Kharchenko et al 



the sky. The system is consid ered as a member of the cluster ( 



fl2005hl in 



1983 



StephensonI ( 



Eggen 



1972, 



Anthony- Twaroa . Il984j ) . The existence of the cluster was suggested by 



I959I ) for the first time. Fur ther photome tric evidence for the 



existence of the cluster was provided by 



Eggen 



(119681 ). However, possible 



members of the cluster, i ncluded BD+36 3317, 



distance modulus values ( lAnthony-Twarogl . Il984( ) and this makes harder to 



lave considerably different 



confirm real members, even the existence of cluster. After all, mean dis- 
tance modulus fo r the c luster is given as 7™. 29 with E{B — V) = 0"'.04 by 



Anthonv-Twarog] (1" 



Kharchenko et al. 



984r) . More recent distance modulus value was given by 



the distance of 373 pc. 



(120051) as 7'".98 (their E 



Kharchenko et al 



(B- V) value is 0"".04, too) with 



(I2OO4J ) calculated the Delta Lyrae 



membership probability of BD+36 3317 as 0.67, 0.91 and 1.0, in terms of 
proper motion, photometry and spatial position, respectively. This calcula- 
tion seems to support the membership of BD+36 3317 to the cluster 
Spectroscopic binary nature of the star has been first noticed by 



Eggen 



(119681 ) via its large radial velocity variation from —9 km s ^ to 17 km s ^. 



Violat-Bordonau and Arranz Her as 



Howev er, eclipses of the system has been discovered by[ 
(I2OO8I ). many years after the discovery of its spectroscopic binary nature. No 
orbital solution via radial velocity measurements has been published up to 



2 



know. At that point, BD +36 3317 has the advantage of to be an echpsing 
binary in terms of estimating its absolute physical properties and distance. 
Furthermore, one can calculate its distance and compare it with the distance 
of the cluster and check the membership of the system to the cluster. For 
those purposes, we obtained standard Johnson U BV observations of the sys- 
tem in 2008 and 2009. By using standard photometric data, we calculated 
E{B — V) value of the system and estimated their spectral type. Although 
the lack of spectroscopic mass ratio and orbital solution is a disadvantage in 
terms of determination of the absolute dimension of the system, we made a 
fair estimation for physical properties of BD +36 3317, including the distance 
of the system, by using the photometric solution. In the next section, we give 
summary of our observations. In section 3 we refine the light elements of the 
system via O — C analysis. In section 4, we give basic photometric properties 
of the system and investigate the effect of interstellar extinction. In section 5 
we present the simultaneous analysis of U BV light curves and our estimation 
for the absolute dimension of the system. In the last section we discuss the 
results. 

2. Observations 

We carried out Johnson U BV observations of the star at Ege University 
Observatory (EUO). We observed the star on 41 separate nights in 2008 
and 2009. Our instrumental setup was 0.3 m Schmidth-Cassegrain telescope 
equipped with uncooled SSP5 photometer. BD +36 3314 was comparison star 
in our observations, while BD +36 3313 was check star. In many observing 
run, when no primary or secondary minimum occurred, we only performed 
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a short observing sequence as S — C2 — CI — WW — S — C2 — CI — 
WW — CI — C2 — S, where 5* is sky, C2 is check star, CI is comparison 
star and V is variable star. For those kind of nights, we used average at- 
mospheric extinction coefficients of EUO in order to correct all differential 
magnitudes in terms oiV — CI, V — C2 and C2 — CI. For other nights, when 
a primary or a secondary minimum occurred, we performed an all night ob- 
serving run. For those nights, we calculate atmospheric extinction coefficients 
via measurements of CI and made all corrections for atmospheric extinction 
on differential magnitudes according to those coefficients. We estimate av- 
erage standard deviations of observations from C2 — CI measurements and 
resulting average standard deviations are 0^.053, O^.OIO and 0^.017 for U, 
B and V filters, resp e ctive 



y. We observed 11 stars from IC4665 cluster 
(iMenzies and Marand . Il996l ) on 29*^ July 2008 and 9 stars from the list of 
19951 ) on I?*'* August 2009, in order to calculate coefficients 
of the transformation of the instrumental system to the standard one. By 
those coefficients, we applied color corrections for all differential measure- 
ments. Then, we directly calculated average standard magnitudes and colors 
of the comparison star from those two nights. Finally, by using standard 
magnitude and colors of the comparison star, we calculate standard magni- 
tude and colors of the variable. 



3. O - C Analysis 

It is not possible to make a comprehensive O — C analysis for BD +36 
3317 since there is not enough minimum time observations. Only the first 



ephemeris of the system is available in literature (jViolat-Bordonau and Arranz Her ad . 
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20081) as 



(HJD) 



Mini 



2,454,437.25921 + r. 30216 x E 



(1) 



where {HJD)Mini is epoch, which corresponds to a time of a primary 
minimum of BD+36 3317, and E is integer echpse cycle number. We have 



aheady had three primary mi nima (Type 



(Type II) in our observations (ISipahi et al. 



) and one secondary minimum 



20091). In Table m we list those 



minima with corresponding O — C values. 



Table 1: The times of light minima of BD+36 3317. In the first column, the errors are 
given for the last digit of the measurements. 



HJD 


E 


0-C 


Filter 


Type 


(24 00000 +) 




(day) 






54652.3522(4) 


50.0 


-0.0025 


UBV 


I 


54667.4148(5) 


53.5 


0.0026 


UBV 


II 


55052.4561(2) 


143.0 


0.0004 


UBV 


I 


55078.2683(4) 


149.0 


-0.0004 


UBV 


I 



Application of linear least squares method to primary and secondary min- 
ima data gives very small amount of corrections in the ephemeris. The new 
light elements and their errors are as follows: 

{HJD)Mini = 2,454,437.2466(30) + 4^^.302162(27) x E . (2) 
For further analysis of light curves, we use those final light elements. 
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4. Basic Photometric Parameters and Interstellar Extinction 

We show phased light and color curves of the system in Figure [H 
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Figure 1: Phase folded and standardized light and color curves of BD+36 3317. 



One can easily notice clear variations in 5 — V in primary and secondary 
minima, which indicate quite large temperature difference among the compo- 
nents of the system. We determined brightnesses and colors of the system in 
maximum light, primary and secondary minima, which would give us some 
hints about components. We list those values in Table HI 

We can inspect effect of interstellar extinction via U BV color - color dia- 
gram, by consi dering the colors at maximu m light. We used UBV standard 



star data from 



Drilling and Landoltl (120001 ) for this purpose. First, we used 
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Table 2: Magnitude and colors of BD+36 3317. 





Max 


Mill I 


Mill II 




(mag) 


(mag) 


(mag) 


V 


8.798 


9.715 


9.053 


U-B 


-0.108 


-0.090 


-0.119 


B-V 


0.066 


0.167 


0.009 



magnitude and colors of maximum light (see Max column in Table H]) by 
assuming E{U — B)/E{B — V) = 0.72 as the slope of the reddening vec- 
tor in UBV color - color diagram. Resulting E{B — V) and values are 
0™.139 and 0"\43, respectively, under the assumption oi R = 3.1. Then, 
de-reddened total color for the system is {B — V)q = -0*^.07. Now we can 
make an estimation for the temperature of the primary component, by as- 
suming that the maximum contribution to the total light comes from the pri- 
mary component. However, during the light curve analysis, we noticed that 
orbital inclination value of the system shows very small variations around 
89.4 degrees which makes the secondary minimum total eclipse. Then, the 
secondary component is completely hidden behind the primary component 
at the middle of the secondary minimum, hence, the magnitude and colors 
of that phase corresponds to the direct measurements of the primary com- 
ponent. Here, we refer the reader to the next section for justification of 
"total eclipse at secondary minimum" case. In this case, we have direct 
measurements of the primary component (third column in Table [2]) and total 
colors and magnitude of the system (first column in Table [2]), therefore we 
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can easily calculate magnitude and colors of the secondary component as V 
= 10^.497, U-B = -0™.053 andB-V = 0^.316. This case is also an ad- 
vantage in terms of more accurate determination of the interstellar reddening 
and intrinsic colors of the components separately. Therefore, we repeated the 
calculation of the interstellar reddening as described above, but this time by 
using the direct measurements of the primary component at the middle of 
the secondary minimum. We list intrinsic colors of the components together 
with the more accurate interstellar reddening and extinction estimation in 
Table H 

Table 3: Intrinsic colors and magnitude of the components of BD +36 3317 together with 
the amount of interstellar reddening and extinction. 





Primary 


Secondary 




(mag) 


(mag) 


Vo 


8.833 


10.277 


{U-B)o 


-0.170 


-0.104 


{B - V)o 


-0.062 


0.245 


E{B - V) 


0.07 


Ay 


0.22 



According to Table |3l colors of the primary component indicates B9 spec- 
tral type while secondary components corresponds about A8 spectral type. 
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5. Analysis of Light Curves 



Under the assumption of total eclipse in secondary minimum, we can 
estimate the effective temperature of the primary component di rectly, which 
is very critical for light curve analysis. We adopted calibration of lGrayl (120051 ) 
for effective temperature estimation of the primary component and resulting 
temperature for {B — V)q = — 0™.062 is Ti = 10750 K with the error of 
(Jti = 470 K. The error of the temperature is calculated from the standard 
deviation of B — V color in our observations. 

Before starting analysis, we chose 0.25 phase as normalisation phase and 
converted all magnitude measurements into normalised flux according to the 



light level at that phase. For 



of th e Wilson - Devinney code ( jWilson and Devinneyl . 



ight curve analysis, we used 2003 version 



1971 



Wilsonl . 



1979, 



19901 ). Photometric properties of the components gives us hints for some 
parameters to reduce the number of free parameters in photometric solu- 
tion. We adopt the gravity brightening coefficients gi = g2 = I and albedos 
Ai = A2 = 1 for stars which have radiative envelopes. We assume syn- 
chronised rotation for both components, so Fi = F2 = 1. We took the 
band-pass dependent (xi 7/19.) and bolome tric {xi^2{bol), yi 2{hol)) limb 



darkening coefficients from 



van Hammd ( 119931 ) by assuming square root law 



( iDiaz-Cordoves and Gimenezl . Il992[ ) which is more appropriate for stars hot- 
ter than 8500 K. Since there is no radial velocity study in literature, we do 
not have any information about mass ratio, which is another critical param- 
eter for light curve analysis. Although it is not an efficient way to determine 
q from photometry in detached systems, we searched for the best solution 
for different q values, starting from q = 0.30 until q = 0.90, by using UBV 
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data simultaneously. Orbital inclination (i), effective temperature of the sec- 
ondary component (T2), Q potentials of primary and secondary (fii, ^^2) and 
luminosity of the primary component (Li) are free parameters in the solu- 
tion. The errors of the best solutions for individual q values are very close 
to each other between q = 0.45 and q = 0.70. At that point, we considered 
mass - luminosity relation as L oc M"^ by using the resulting absolute lu- 
minosities at the end of the solution (see later results in this section). We 
repeated solutions for many q values between q = 0.45 and q = 0.70. In most 
cases, L oc relation indicates the q value close to 0.65, hence, we finally 
accepted q = 0.65 and applied a final light curve solution. In Table HI we 
give final light curve analysis results. We note that the error of T2 is internal 
to the Wilson - Devinney code and its error should be similar to the error of 
T2. In the table, (^1^2) denotes average of three fractional radius values {pole, 
side and back values in solution output) relative to the semi-major axis of 
the orbit, for corresponding component. We give the normalised fluxes for 
U, B, and V filters and corresponding theoretical light curves for our final 
solution in Figure |2j 

In Figure |3l we zoom to the primary (left panels) and secondary (right 
panels) minima to show the shapes of the eclipses. One can notice that the 
secondary minimum is certainly total eclipse which lasts for a short phase 
range. At the primary minimum, we observe non-fiat bottomed light varia- 
tion which shows the effects of annular eclipse and limb darkening together. 

Comparison among photom etric solution results , mass - luminosity rela- 



tion and evolutionary models of 



Girardi et al. 



(|2000[ ) indicates that the mass 
of the primary Mi is close to the 2.5 Mq, therefore we assume Mi = 2.5 Mq. 
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Table 4: Light curve analysis results of BD +36 3317. Errors are given in parenthesis. 



M2/M1) 


0.65 (fixed) 


Ti{K) 


10750 (fixed) 


91 = 92 


1.0 


Ai = A2 


1.0 


Fi — F2 


1.0 


i (°) 


89.61(11) 


T2{K) 


7711 (10) 




10.571 (25) 




8.997 (23) 


Li/{Li+L2)u 


0.851(17) 


Li/(Li+L2)ij 


0.831(15) 


Li/{Li+L2)v 


0.789(15) 


Xi{hol), X2{hol) 


0.558, 0.215 


yi{bol),y2{bol) 


0.172, 0.525 




0.082 , 0.590 




-0.058 , 0.846 




-0.047 , 0.723 


X2,y2{U) 


0.184 , 0.646 


X2,y2{B) 


0.105 , 0.822 


X2,y2(V) 


0.096 , 0.711 


(ri),(r2) 


0.1009(3), 0.0832(2) 


rms 


0.011 
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Figure 2: Representation of observational data (points) in terms of normalised flux and 
theoretical solution (continuous line). 

This assumption makes the mass of the secondary component M2 = 1.6 Mq 
according to the q value. 

Since we estimated the masses of the components, we can go one step 
further and calculate the absolute parameters and the distance of the sys- 
tem. By the aid of Kepler's third law, we can calculate semi-major axis of 
the orbit a. After that point, one can easily calculate absolute radii of the 
components by using average fractional radii in Table |H Now we have effec- 
tive temperatures and absolute radii of the components and we can calculate 
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Figure 3: Close look to the primary and secondary minima in each filter. In secondary 
minimum, one can notice the short total eclipse. 



their luminosities in solar unit via Stefan — Boltzmann law, by using Tq = 
5770 K. We list absolute parameters of the system in Tabled Finally, we can 
calculate the distance of the system by using primary and secondary compo- 
nent separately, via their photometric and absolut e properties. In distance 
calculation, we adopted bolometric corrections from I Gray! (120051 ) . Photomet- 
ric properties of the primary component leads to a distance of 353 pc. We 
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only use primary component in order to estimate the distance since its signal 
is very strong relative to the secondary component. 



Table 5: Absolute Physical Properties of BD +36 3317. Solar Mm value of 4.^74 is used 
to calculate Mhoi values of the components. 



Parameter 


Primary 


Secondary 


Mass (Mq) 


2.5 


1.6 


Radius (Rq) 


1.8 


1.5 


Luminosity (L©) 


39 


7 


Mbol (mag) 


0.76 


2.62 


log{g) (cgs) 


4.32 


4.31 



6. Summary and conclusions 

We presented standard Johnson UBV photometry of the Algol type 
eclipsing binary BD+36 3317 with a fairly good phase coverage and rea- 
sonably accurate observational data. During observations, we obtained three 
primary and one secondary minima. We refined the epoch and the period of 
the system by applying hnear least squares method to the timings of those 
light minima. We determined standard colors and magnitude of the system 
in maximum light, primary minimum and secondary minimum. During the 
analysis, we noticed that the secondary minimum is total eclipse, which is an 
advantage in analysis and means that the measurement at that phase cor- 
responds to direct measurements of the primary component. Using this ad- 
vantage, we first determined interstellar reddening and extinction via direct 



14 



measurements of the primary component. Then, we were able to calculate 
de-reddened colors and magnitudes of the components, separately. This case 
was an another advantage in order to make a more accurate estimation of Ti, 
hence enabled us to reduce the number of free parameters more reliably in 
simultaneous UBV light curve solution. Photometric solution justified that 
the secondary minimum was total eclipse. 

Lack of an orbital solution based on radial velocity measurements pre- 
vented us from determining the accurate absolute parameters of the system 
and their uncertainties. Hence, we can only make a rough estimation for the 
uncertainties and check how those uncertainties affect our results. 

If we assume that the q is between 0.6 and 0.7 and Mi is between 2.4 Mq 
and 2.6 Mq, then we can calculate a lower and upper limit for a via Kepler's 
third law. Those lower and upper limits of a put -Ri between 1.76 Rq - 1.84 
Rq when we consider the average fractional radius of the primary component 
from photometric solution. Same calculation puts R2 between 1.45 Rq - 1.52 
Rq. a similar method can be used for the luminosities of the components 
by using Stefan — Boltzmann law and Tq = 5770 K. Assuming a lower and 
upper limits for Ti and T2 via ctt^ = 470 K, we can calculate the ranges of 
Li and L2 as 31 Lq - 49 Lq and 5 - 9 Lq, respectively. 

We show preliminary plots of the components in log Tg// - log L plane in 



Figure [H Evo 



Girardi et al. 



utionary tracks for solar abundance {Z = 0.019) comes from 



( I2OOOI ). The components of the system seem close to the ZAMS 



and in good agreement with solar metal abundance. However, spectroscopic 
analysis is necessary to revise or refine it. 

We can make an estimation for the uncertainty of the distance in a sim- 
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Figure 4: Positions of the components (fille d stars) in log T^/f - log L plane. All tracks 



ilar way, as described above. If we assume the same ranges for q, Mi and 
Ti as in previous uncertainty estimations, we can calculate the range of ab- 
solute bolometric magnitude of the primary component as 1'".01 - 0™.52. 



Here, we adopt bolometric corrections from iGrayl ( l2005l ) for corresponding 
temperature limits in order to calculate limit visual absolute magnitudes 
(My) of the primary component. Using de-reddened V magnitude and My 
limits of the primary component in distance modulus, we can calculate the 
range of the distance as 329 - 377 pc and this leads a mean value of 353 pc 
which is our estimation in previous section. The uncertainty in the distance 
might be slightly exaggerated since we make a rough uncertainty estimation 
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for related parameters. However, our distance estimation differ s only by 20 



2005|). This con- 



pc from distance value of the cluster flKharchenko et al. . 
firms the membership of BD+36 3317 to Delta Lyrae cluster. Nevertheless, 
comprehensive spectroscopic study of the system, in terms of radial velocity 
measurements, would help to refine or revise the physical properties and dis- 
tance of the system. Further spectroscopy would also give a chance to check 
the metal abundance of the system which also contains some hints about the 
nature of the cluster. 
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